Persistent localized states for a chaotically mixed bistable reaction 
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We describe the evolution of a bistable chemical reaction in a closed two-dimensional chaotic 
laminar flow, from a localized initial disturbance. When the fluid mixing is sufficiently slow, the 
disturbance may spread and eventually occupy the entire fluid domain. By contrast, rapid mixing 
tends to dilute the initial state and so extinguish the disturbance. Such a dichotomy is well known. 
However, we report here a hitherto apparently unremarked intermediate case, a persistent highly 
localized disturbance. Such a localized state arises when the Damkohler number is great enough to 
sustain a "hotspot" , but not so great as to lead to global spread. We show that such a disturbance is 
located in the neighborhood of an unstable periodic orbit of the flow, and we describe some limited 
aspects of its behavior using a reduced, lamellar model. 

PACS numbers: 82.20.-w, 82.40.Bj 



I. INTRODUCTION 

The progress of chemical or biological processes mixed 
by fluid flows is of great environmental and industrial 
significance, in a wide variety of contexts (see jjj. L3- Lj- L| 
for some broad applications; in this paper we shall, for 
brevity, mostly use the language of chemical processes 
in our discussions). The underlying processes often pos- 
sess multiple possible long-time states, characterized by 
whether some reaction proceeds or dies out. The ulti- 
mate fate of a localized initial disturbance can depend 
subtly on the combined influences of advection by the 
fluid, diffusion of the reactants, and the reaction kinet- 
ics. 

This paper concerns the influence of chaotic, laminar 
fluid mixing on the course of a bistable chemical reac- 
tion in a closed flow domain, fn a spatially homogeneous 
environment, such a reaction tends to one of two steady 
states (which we describe below as the "quenched" and 
"excited" states) according to whether or not the initial 
reactant concentrations exceed some threshold. However, 
the large-time behavior of the reaction from a spatially 
localized disturbance subject to fluid mixing is a more 
subtle question, since the mixing has two competing ef- 
fects on the course of the reaction. On the one hand, the 
mixing promotes the chemical reaction, since, according 
to the kinematics of the fluid motion, interfaces are ex- 
ponentially stretched. However, the mixing also tends to 
suppress the reaction, since the repeated stretching and 
folding that is characteristic of chaotic flows generates a 
complicated filamentary pattern; if the thinning of the 
filaments is sufficiently rapid, it may ultimately quench 
the reaction, through diffusive dilution of the reactants. 
The loss of a stable excited state as the mixing grows 
in strength (or as the reaction rate is reduced) may be 



viewed as a saddle-node bifurcation, and has been identi- 
fied in a number of reaction schemes 0, @, 0, U E3> EI • 
For example, Neufeld et al. note that for both bistable 
and autocatalytic reactions there are two regimes: when 
mixing is slow, localized excited perturbations propa- 
gate throughout a closed domain; by contrast, under 
rapid mixing, localized perturbations decay to the uni- 
form quenched state. In the propagative case, perturba- 
tions spread in the form of filaments. For an excitable 
medium, Neufeld et al. [jj describe three scenarios for a 
closed flow: in the first, there is eventually a global ex- 
citation, which is spatially incoherent, so different parts 
of the domain are excited at different times; in the sec- 
ond, there is again a global excitation, but this time it 
is spatially coherent, and ultimately the system homo- 
geneously decays to the unexcited state; in the third, 
the excitation enjoys no significant propagation before 
its ultimate decay. Kiss et al. 01 investigate a chem- 
ical reaction scheme capable of supporting both steady 
and oscillatory homogeneous states, subject to chaotic 
mixing. They find that ultimately the system reaches 
a spatially homogeneous state, which may be of either 
steady or oscillatory type, depending on the parameters. 

We shall show here the existence of a third regime — 
beyond the global quenched and excited states — com- 
prising a persistent localized disturbance. For definite- 
ness, we consider the model bistable reaction scheme in 
which a chemical concentration A(x,t) satisfies, in di- 
mcnsionless form, 



dA 
~dt 
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Pe" 1 V 2 A + Da^(l-yl)( J 4-a), (1) 
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where < a < 1/2 [E 0, 13. Here Pe = UL/D is the 
Peclet number, Da = kC 2 /UL is the Damkohler number 
and u(x,t) is the dimensionless fluid velocity (here U, 
L and C are scales for the fluid velocity, lengths and 
chemical concentrations, respectively; D is the diffusion 
coefficient and k is the reaction constant). In view of the 
need to carry out well resolved simulations at moderately 
high Peclet number, we shall consider a two-dimensional 
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problem, in which it is the so-called "sine flow" . This 
choice of flow permits a particularly efficient and accurate 
spectral numerical solution of the problem [THl ITfij . 

Let us now discuss the behavior of (JJJ. First consider 
the reaction in the spatially homogeneous case. Then if 
the initial value of A exceeds the threshold a, A — > 1 
at large times. By contrast, if A is initially less than 
the threshold, then A — > at large times. We expect a 
similar scenario to hold in if the initial state is non- 
homogeneous, but the mixing is sufficiently rapid: an 
initial state is first rapidly averaged in space, then the 
long-time behavior of the system depends on whether 
or not, in this averaged (nearly spatially homogeneous) 
state, A exceeds a. In the opposite limit of no fluid mo- 
tion, an initial localized disturbance spreads or becomes 
extinct depending on the initial state, but if the distur- 
bance takes the form A — Aq (a constant) in some finite 
region f2, and A = outside O, then the condition for the 
disturbance to spread is roughly that Aq should exceed a. 
We expect a similar scenario for slow fluid mi xing . These 
considerations, and the results of Refs . ml l7L I9L suggest 
that, in a closed flow, the fate of a localized disturbance 
in 0) is as follows: the system tends to a uniform state 
at large time, with A = (for sufficiently rapid mixing 
or slow reaction or weak initial perturbation) or A = 1 
(for sufficiently slow mixing or rapid reaction). Further- 
more, given an initial disturbance that leads to the lat- 
ter globally excited state for given values of Da and Pe, 
one would expect that a sufficient reduction in Da would 
cause the reaction to proceed instead to the quenched 
state. However, we show here that there is an additional, 
intermediate scenario, in which the perturbation appears 
to remain localized for all time, spread along a segment 
of the unstable manifold of an unstable periodic orbit 
of the flow. We investigate this localized state through 
numerical simulation of in particular exploring its 
dependence on Damkohler and Peclet numbers, and on 
the initial state. Finally, we shed some limited analytical 
light by consideration of a reduced, lamellar model that 
represents the balance between spreading of the hotspot 
due to reaction and diffusion, and compressive effects of 
the flow. 



II. SIMULATIONS 

In this section we describe some numerical simulations 
of (JTJ. Stirring of the fluid is accomplished by the "sine 
flow", which is time-periodic with period T; then |14| 




(sin 2th/, 0) for mT <t< (m + \)T, 
(0,sin27ra;) for (m + |)T < t < (to + 1)T, 



for to = 0, 1, 2, . . .. The simulations take place in a dou- 
bly periodic square box < x,y < 1. We consider the 
"globally chaotic" case T — 1.6, for which the Poincare 
section appears to the eye to have no significant regular 
islands. The initial condition is given by the Gaussian 




FIG. 1: (a) Time evolution of the spatially averaged concen- 
tration (A) for Pe = 10 4 , from the initial condition with 
Q — 1, for various Damkohler numbers, as indicated. 

concentration profile 

A(x, y, 0) = exp [-100tt((z - 1/4) 2 + (y - l/4) 2 )/Q] , 

(3) 

for < x, y < 1, so that the initial spatially averaged 
concentration is (A) = 0.01 x Q at t = 0. In most sim- 
ulations, Q = 1, although we also briefly consider other 
values of Q, in the range < Q < 2. The choice to cen- 
ter the initial disturbance on the point (0.25, 0.25) was 
arbitrary, but, as we shall see, it proves serendipitous. 

We simulate Q subject to the flow field J2J and the 
initial state © using a pseudospectral code, which is de- 
scribed in Ref. ^] (see also Ref . fLsT) . The Peclet number 
in our simulations is limited to 10 5 , to ensure well re- 
solved results with our available computational resources. 
Most simulations up to Pe — 10 4 arc carried out using 
512 Fourier modes in each of the x and y directions, and 
a time step of 2 x 10~ 3 ; unlike the simulations reported 
elsewhere |l6j . which exclusively use the second-order 
time-stepping scheme ETDRK2, we mostly use here the 
fourth-order scheme ETDRK4 (see Ref. for details of 
these two schemes). 

In all the numerical simulations below, the sine flow is 
fixed to have period T = 1.6, in order to avoid the appear- 
ance and disappearance of regular islands with varying T. 
Thus our investigation of the effects of fluid mixing is car- 
ried out by varying the Peclet number (smaller and larger 
values of this parameter corresponding, respectively, to 
more rapid and slower mixing). We also investigate the 
influence of the reaction rate and the initial amount of 
A upon the long-time behavior of the system through 
varying the Damkohler number Da and Q, respectively. 



A. Numerical results 

In all our simulations, we set a = 0.25. Thus for suffi- 
ciently rapid mixing the reaction is quenched (since the 
initial averaged concentration (A) = 0.01Q < 0.25), and 
for slow enough mixing the disturbance spreads through- 
out the domain (since the peak initial concentration 
A max = 1 > 0.25). 
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FIG. 2: The circles are centered on a period-6 periodic orbit 
of the sine flow @, plotted at times T, 2T, 6T (recall 
that the line x = is identified with the opposite edge x — 1) . 
Also shown are the contours A(x,y, t) = 0.5 at corresponding 
times, for a hotspot found at Da = 6 and Pe = 10 , indicating 
that the hotspot is located in the vicinity of the period-6 orbit. 



We first set Pe = 10 4 and Q = 1 and consider the 
effect of varying Da. Fig. ^ shows some illustrative re- 
sults. For Da = 6, the reaction is too slow to sustain 
itself and A — > everywhere as t — > oo. For Da = 12, 
the reaction is fast enough to overcome the diluting influ- 
ence of the fluid mixing; in this case A — v 1 everywhere 
as t — > oo. The choice Da = 8 provides an intermediate 
case, for which the reaction proceeds to neither homoge- 
neous state as t — > oo, but rather generates a persistent 
time-dependent state in which A = in most of the flow 
domain, but A « 1 in a small time-dependent region. 
This region turns out, upon further examination, to lie 
in the vicinity of a period-6 orbit of the underlying flow 
field u. The localized state appears to be time-periodic, 
with a period of 6T, while the spatial mean concentra- 
tion (A) appears to have period 3T. We have tracked 
this solution for many tens of periods, and it is robust. 
A further simulation, at Da = 9, at first appears to pro- 
duce another (somewhat larger) localized "hotspot" , but 
the hotspot does not approach a time-periodic state; in- 
stead it gradually grows, eventually filling the entire flow 
domain, so that (A) — > 1. 

The localization on an unstable periodic orbit is read- 
ily understood in qualitative terms, as follows. In the 
vicinity of such an orbit there are directions of local 
compression and stretching. In the compressive direc- 
tion, a balance is maintained between the spreading of 
the reaction-diffusion front and the compressive veloc- 
ity field H Q, fl3| . resulting in a finite "width" for the 
hotspot. The hotspot is also stretched along the unsta- 
ble manifold of the periodic orbit; growth in its length 
seems to be moderated by the inability of the reaction to 
sustain a thin filamentary structure in the vigorous fluid 
mixing. 

We have located this unstable period-6 orbit, us- 
ing the stabilization algorithm of Pingcl, Schmelcher 
and Diakonos 01 ( see Fig. [SJ: at times mT, where 
m is an integer, the points {x m ,y m ) on the orbit are 
••• -> (0,0.2935) -> (0.7032,0.5) -> (0.7032,0.7065) -> 
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FIG. 3: Long-time behavior of Q, subject to |JSJ and with 
Q = 1, in Da-Pe parameter space. The three long-time fates 
are: (A) — <> (*); a localized hotspot (4-); (A) — > 1 (x). 



(0,0.7065) -> (0.2968,0.5) -> (0.2968,0.2935) 

We note that, allowing for the spatial periodicity of 

the flow in both x and y, we have {x m +z,ym+?,) 

1 x m , 1 — y m )i corresponding to a rotation of 7r radi- 
ans about the center of the flow domain. A corresponding 
symmetry may be observed in the concentration contours 
in Fig. [3 (following from the point symmetry of the ve- 
locity field u about the center of the flow domain), ex- 
plaining why (A) has period 3T rather than 6T. Recall 
that our initial condition is centered on the point 
(0.25,0.25), which is close to the point (0.2968,0.2935) 
on the period-6 orbit; this fortunate proximity underlies 
our finding the hotspot. 

We compute that the Jacobian matrix associated with 
the period-6 orbit has eigenvalues 18.73 and 0.053. In the 
next section we shall use this information regarding the 
stretch rate close to the periodic orbit in the development 
of a reduced model for the hotspot. 

Fig. |3 summarizes the results of a large set of sim- 
ulations, all at Q — 1, for different Damkohler and 
Peclet numbers. Since the goal of the simulations was 
to determine the boundaries in parameter space between 
the three different long-time behaviors, we have not at- 
tempted to fill in data where the behavior is (presum- 
ably) predictable. (It should be noted that each simula- 
tion currently takes on the order of a day or two to run 
on a desktop PC.) Gaps in the data towards the lower 
right-hand part of the figure reflect the increased length 
of runs required there for the long-time behavior to be- 
come apparent. 

Consider first the parameter dependence to the left of 
Fig. [3J if we start in the region of hotspots and decrease 
either Da or Pe then we move into the (A) — > region. 
This is because either the reaction rate becomes too low 
to sustain the hotspot, or the initial rate of spreading 
of the chemical becomes too great for an adequate level 
of material to be delivered to the period-6 orbit in the 
first instance. Hence the system is quenched. To the 
right of the figure, if we start in the hotspot region and 
increase Da then we move into the (A) — > 1 region, be- 
cause the reaction rate becomes too great for the hotspot 
to remain localized. Similarly, if we instead reduce Pe, 
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FIG. 4: Time evolution of (A) for Da = 8 and Pe = 10 4 , 
from the initial condition for various values of Q. (The 
value of Q corresponding to each curve may be determined 
as 100 times the initial value of {A}.) For Q — 2,1.5,1.4, 
A — * 1 everywhere at long time; for Q = 0.38, 0.35, 0.3, ^4^0 
everywhere at long time. Intermediate values of Q lead to the 
(same) hotspot solution; such results are shown for Q = 0.39 
and for Q from 0.4 to 1.3 in steps of 0.1. 



then again the hotspot cannot remain localized, due to 
its enhanced rate of spreading through diffusion, and the 
system is globally excited. The parameter dependence of 
the system towards the lower right-hand side of the data 
in Fig. |3| is rather subtle, since a reduction in Pe may 
lead either to the quenched state or to the excited state. 
Indeed, we find that our numerical simulations in this re- 
gion of parameter space often involve a transient, slowly 
evolving hotspot before eventually selecting the ultimate 
quenched or excited state. 

We next discuss the influence of the initial amount of 
A upon the long-time behavior of the system. To this 
end, we have carried out simulations with Pe = 10 4 and 
Da = 8 for a range of values of Q. As Fig. ^ shows, 
Q > 1.4 leads to the globally excited state, Q < 0.38 
leads to the globally quenched state, and intermediate 
values of Q lead to a hotspot at long times. Although 
it is not apparent from the figure, which focuses on rela- 
tively short-time results, the long-time hotspot solution 
is the same for all of these simulations with Q in the 
appropriate range (this is evident in a sample of longer 
simulations that we have carried out, but which are not 
reported here). In other words, for these simulations the 
ultimate hotspot solution depends only on the system 
parameters (i.e., T, Pe and Da); the role of the initial 
condition is merely to determine which of the three long- 
time states is selected. 

We now turn to the Peclet number dependence of the 
hotspot, characterizing the rate of fluid mixing. In Fig. [5] 
we fix the reaction rate Da = 6 and Q — 1 and vary Pe. 
Here we find that at Pe = 10 4 , as we have described 
above, the relatively rapid mixing quenches the reac- 
tion, so that A — > everywhere. At Pe = 10 5 /3 and 
at Pe = 10 5 , by contrast, the mixing is slower, and we 
again find a persistent hotspot. If we accept the qualita- 
tive picture above, that the width of the hotspot is deter- 
mined by a balance between compression and reaction- 



FIG. 5: Time evolution of (A) for Da = 6, from the initial 
condition 0, for various Peclet numbers, as indicated. 

diffusion, while its length is determined by mixing in the 

I/O 

bulk then the width should scale as Pe ' for large Pe, 
and the length should remain roughly independent of Pe. 
Then the mean concentration (^4) should scale as Pe -1 ^ 2 . 
For the two hotspots indicated in Fig. [5J this scaling is 
borne out, although admittedly over a very short range: 
we find, for the minimum value of (A) over the cycle, 
that Pe 1/2 (A) min = 0.98 and 1.02 at Pe = 10 5 /3 and 
10 5 , respectively. (These values would be equal if the 
Peclet number scaling of (A) were exact.) The corre- 
sponding values of the maximum mean concentrations 
over the cycle are Pe 1 ^ 2 (A) max = 1-15 and 1.10. If the 
qualitative argument holds, then one should expect that 
as Pe is increased beyond 10 5 , the hotspot persists, with 
(A) continuing to scale with Pe -1 / 2 . Unfortunately we 
are unable to simulate reliably beyond Pe = 10 5 , so we 
are unable to confirm this prediction. 

III. LAMELLAR MODEL 

There has recently |lll[lllll[]3[lJ[ll[]llI3 

been significant progress in understanding the dynamics 
of reaction-advection-diffusion problems such as Q by 
analysing reduced "filamental" or "lamellar" models [21| , 
which concern the simpler problem of the evolution of the 
reaction in a one-dimensional setting, in our case 

8 A 8 A 8 2 A 

^-^=Pe- 1 ^ J +DaA(l-A)(A- a ). (4) 

Here £ represents a coordinate in the compressive direc- 
tion near an unstable periodic orbit of the sine flow. This 
is the "Lagrangian filament slice model" , which has been 
thoroughly reviewed by Tel et al. [22J, largely, but not 
entirely, in the context of open rather than closed flows. 

This model might be expected to provide at least a 
partial picture of the behavior of the hotspot. To see 
this, we first note that the hotspot is a thin structure, 
which appears, from the two-dimensional numerical sim- 
ulations of Sec. Ill Al to be approximately concentrated 
around the unstable manifold of an unstable periodic or- 
bit. Thus we approximate the flow in the neighborhood 



5 



0.8 - 
0.6 - 




-0.02 -0.01 0.01 0.02 

FIG. 6: Cross-section through the hotspot at Pe = 10 and 
Da = 6. Dashed line: from Q, with £ being a coordinate 
roughly normal to the hotspot. Solid line: corresponding pre- 
diction from Q|. 



of the orbit by its linearization. The key balance in de- 
termining the existence or otherwise of the hotspot is 
that between the compressive flow in the direction of the 
stable manifold (represented by the term — \£A^ in Q) 
and the tendency of the hotspot to spread as a reaction- 
diffusion front (represented by the right-hand side of 
and so we reduce consideration to one-dimensional vari- 
ations along the stable manifold, and ignore variations 
along the unstable manifold. This remains a complicated 
model, and so we take the further gross liberty of assum- 
ing that the flow field along the stable manifold is time- 
independent, with a uniform and constant rate of com- 
pression A. For example, for the period-6 orbit, the value 
of A is chosen such that exp(— 6AT) corresponds to the 
compressive eigenvalue of the Jacobian matrix of the pe- 
riodic orbit (i.e., A = (- log 0.0534)/ (6T) = 0.3052). We 
are thus led to as a simplified one-dimensional model. 
Since nontrivial solutions to are localized around the 
origin, we may make the further simplification that the 
corresponding spatial domain is — oo < £ < oo. 

Let us now briefly recall some relevant known results 
for (0J. An exact characterization of the behavior of a 
single reaction-diffusion front in (0J is possible in the 
absence of the advection term — A£^4j [(J |2^, [2j|: such 
a front between A = 1 as ^ -> - oo and A = as 
£ — > oo moves with velocity v = (Da/2Pe) 1 / 2 (l — 2a). 
Furthermore, its profile is readily determined to be A = 
{1 + cxp[/i(£ — vt)]} 1 (up to some shift in time or space 
origin), where /x = (DaPe/2) 1 / 2 . 

In the presence of the advection term, such an exact 
characterization of solutions to (0} is no longer possi- 
ble, although (large-Da) asymptotic and numerical treat- 
ments are available |13| . Of particular interest in the 
present context is the behavior of solutions to on 
— oo < £ < oo from localized initial disturbances, which 
is well documented H IH El El: if Da/A is less than 
some threshold f(a) then all initial states decay to A = 0; 
however, if Da/A > f(a) then an above-threshold initial 
state generally tends to a nontrivial, stable, steady lo- 
calized state as t — > oo. The transition at Da/A = f(a) 
is a saddle-node bifurcation For a = 0.25, as 



in our two-dimensional simulations, we find this saddle- 
node bifurcation to be at /(0.25) = 11.13 (by solving Q 
numerically and tracking the stable steady state as Da 
is reduced). Since A = 0.3052, we thus conclude that 
no hotspots should be possible in the two-dimensional 
simulations for Da < 3.4. It is certainly the case that 
we have observed no hotspots for such small Damkbhler 
numbers; indeed, the lowest Damkohler number for which 
we have observed hotspots is 6. Given the small sample 
of parameter space provided by our simulations and the 
gross assumptions made in arguing for this seems to 
represent reasonable agreement between this model and 
the full simulations. Further quantitative comparisons 
may be made by examining the structure of the hotspot 
itself. To this end, Fig. shows a cross-section of the 
concentration through the hotspot, from two-dimensional 
simulations of and from the model (@J. The gen- 
eral concentration profile is captured well by the model, 
although the hotspot itself is rather narrower than the 
one-dimensional model prediction. The one-dimensional 
model also predicts that only sufficiently wide Gaussian 
initial conditions will generate a hotspot, for supercriti- 
cal Damkohler numbers |l3[ . This dependence upon the 
initial condition is qualitatively consistent with the full 
system, as illustrated in Fig. 21 

Of course, the one-dimensional model has only lim- 
ited predictive power for the full two-dimensional system. 
While the saddle-node bifurcation of the model seems 
to capture the demise of the hotspot as the Damkbhler 
number is decreased, as illustrated in Fig. ^ there is no 
mechanism in the model for global excitation (instead the 
model predicts a hotspot whose width grows like Da ' 
for large Da [L|)- In other words, the model predicts be- 
havior qualitatively consistent with the left-hand bound- 
ary of the hotspot region in Fig. El but not the right-hand 
boundary. A further shortcoming of the model is that 
its dependence upon the Peclet number can be removed 
by a straightforward scaling of £, unlike the full two- 
dimensional problem. As a consequence, the saddle-node 
bifurcation in is independent of Pe and, for a given 
Damkohler number, the stable steady-state solution sim- 
ply becomes thinner as the Peclet number is increased, 

1/2 

with width proportional to Pe ' . Such behavior is in 
contrast to the precipitate loss of the hotspot illustrated 
in Fig.|5]for Pe = 10 4 , and the left-hand boundary of the 
hotspot region in Fig. [3J 



IV. CONCLUSIONS AND DISCUSSION 

For a bistable chemical reaction mixed in a chaotic 
laminar fluid flow, we have demonstrated numerically lo- 
calized states that seem to persist indefinitely. In a sweep 
of parameter space, these localized states provide a case 
intermediate between global quenching and global excita- 
tion of the reaction, as the Damkohler number is varied. 
They are associated with unstable an periodic orbit, and 
their existence may partially be understood in terms of 
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a simple lamellar model. 

The long-time behavior of the system is strongly de- 
pendent on the size of the initial disturbance: too large or 
small and the system approaches, respectively, globally 
excited or quenched states; only for initial disturbances 
of intermediate size can a hotspot be found. 

The broad influence of the Damkohler number upon 
the long-time state of the system is straightforward to 
summarize: if we begin with parameter values that lead 
to a hotspot, then a sufficient increase or decrease in 
the Damkohler number takes the system instead to the 
globally excited or quenched states. 

The influence of the Peclet number is slightly more 
complicated. Again suppose that we start with parame- 
ter values leading to a hotspot, and consider the effects 
of reducing the Peclet number (i.e., enhancing diffusion). 
For lower reaction rates (i.e., smaller Damkohler num- 
bers), we find that the enhanced initial smearing of a 
localized initial disturbance results ultimately in a glob- 
ally quenched state. When the reaction rate is greater, 
however, the nascent hotspot is able to sustain itself even 
in the presence of the enhanced diffusive spreading, but 
is unable to remain localized when the Peclet number 
is reduced sufficiently; thus the system approaches the 
globally excited state. 

Since we are unable accurately to simulate Peclet num- 
bers greater than approximately 10 5 , any large-Peclet 
number trends are inaccessible to our numerics. Nev- 
ertheless, we conjecture that a hotspot generally persists 
as the Peclet number is increased (all other parameters 
being fixed), with the amount of chemical present seal- 

I/O 

ing as Pe ' . Of course, an exception to this suggestion 
would arise if one chose an initial condition for which 
diffusion was essential in transporting enough of an ini- 
tial patch of chemical to the vicinity of the periodic orbit 
upon which the hotspot would ultimately be attached. 
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